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A total of 137 upward stopping muons of minimum energy 1.6 GeV are ob- 
served by Super-Kamiokande during 516 detector live days. The measured muon 
flux is 0.39 ± 0.04(stat.) ± 0.02(syst.) x 10^"cm"^s"^sr~^ compared to an expected 
flux of 0.73 ± 0.16(theo.) x 10^^'^cm^'^s^^sr^^. Using our previously-published mea- 
surement of the upward through-going muon flux, we calculate the stopping/through- 
going flux ratio TZ, which has less theoretical uncertainty. The measured value of TZ 
=0.22±0.02(stat.) ± 0.01(syst.) is significantly smaller than the value 0.37t° °4(theo.) 
expected using the best theoretical information (the probability that the measured TZ 
is a statistical fluctuation below the expected value is 0.39%). A simultaneous fitting 
to zenith angle distributions of upward stopping and through-going muons gives a re- 
sult which is consistent with the hypothesis of neutrino oscillations with the parameters 
sin^2^ > 0.7 and 1.5xl0"^<Am^<1.5xl0"^ eV^ at 90% confidence level, providing a 
confirmation of the observation of neutrino oscillations by Super-Kamiokande using the 
contained atmospheric neutrino events. 



PACS numbers: 14.60.Pq, 96.40.Tv 



keywords: neutrino oscillations, Super-Kamiokande, upward stopping muons, upward throughgo- 
ing muons, atmospheric muon neutrinos, cosmic rays 
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Energetic atmospheric or i/^ interact with the rock surrounding the Super-Kamiokande ( "Super-K" ) 
detector and produce muons via weak interactions. For downward-going particles, the rock overburden is 
insufficient to prevent cosmic-ray muons from overwhelming any neutrino induced muons, but upward- 
going muons are or induced because the entire thickness of the earth shields the detector. Muons 
energetic enough to cross the entire detector are defined as "upward through-going muons" , and have 
been discussed in a previous article The typical energy of the parent neutrinos is approximately 
too GeV. Upward-going muons that stop in the detector are defined as "upward stopping muons" , and 
come from parent atmospheric neutrinos with a typical energy of about 10 GeV. These energy spectra are 
shown in Fig. 0. Neutrinos arriving vertically from beneath the detector travel roughly 13,000 km from 
their point of production, while those coming from near the horizon originate only ^--^500 km away. Thus, 
observation of upward-going muons provides a relatively pure sample of muon neutrinos, with a wide 
range of path lengths, allowing tests of possible i/^ disappearance due to flavor neutrino oscillations 

For atmospheric neutrinos with interaction vertices inside the detector's fiducial volume (referred to 
as "contained" events), Super-K measures a low I'^/i^e ratio |^ and clearly observes a strong zenith 
angle dependence of the i^f^ events [Q . This muon neutrino disappearance |^] is consistent with ^ i^r 
oscillations, while the up-down symmetry seen in the Ve events rules out significant Ve appearance, in 
agreement with recent results from the CHOOZ experiment The z^^ ^ i^t oscillation hypothesis 
suggested to explain the Super-K contained event results is also consistent with anomalous upward 
through-going muon zenith angle distributions observed by Kamiokande 0, MACRO Q and Super- 
Kamiokande Q. The upward stopping muons are the remaining class of atmospheric neutrino events 
with which Super-K can test the oscillation hypothesis. It should be noted here that the parent neutrino 
energies of upward stopping muons are similar to those of multi-GeV fully-contained (FC) and partially- 
contained (PC) events discussed in Thus, a similar deficit of upward stopping muons to that of 
upward-going multi-GeV FC and PC events is expected if the oscillation hypothesis is the cause of the 
multi-GeV FC and PC event deficit. The main difference between the multi-GeV FC and PC event 
sample and the upward stopping muon sample is that roughly 80% of the upward stopping muons are 
generated in the surrounding rock, which has different neutrino interaction cross-sections and systematic 
uncertainties than does water. 

The Super-K detector is a 50 kton cylindrical water Cherenkov detector. To reduce the cosmic-ray 
muon background, the detector was constructed ~1000 m (2700 m.w.e.) underground at the Kamioka 
Observatory, Institute for Cosmic Ray Research, the University of Tokyo, in the Kamioka-Mozumi mine, 
Japan. The detector is divided by an optical barrier instrumented with photomultiplier tubes ( "PMT" s) 
into a cylindrical primary detector region (the Inner Detector, or "ID") and a surrounding shell of water 
(the Outer Detector, or "OD" ) serving as a cosmic-ray veto counter. Details of the detector and general 
data reduction procedures can be found in reference . The data used in this analysis were taken from 
Apr. 1996 to Jan. 1998, corresponding to 516 days of detector livetime. 

An upward-going muon is defined as a track that appears to enter the detector from the rock, and 
reconstructs as traveling in the upward direction. Thus, PMT activity in the OD at the muon's entrance 
point is required. The total cosmic-ray muon rate at Super-K is 2.2 Hz, of which a few percent are 
stopping muons, and the great majority are downward-going. The trigger efficiency for a muon entering 
the ID with momentum more than 200 MeV/c is ■~100% for all zenith angles. 

Muons which leave entrance signal clusters in the OD with no corresponding exit cluster are regarded as 
stopping. A neutrino interaction inside the water of the OD also produces such a signal, so some fraction 
(~ 20%) of the upward-going stopping muon sample is estimated to originate in the OD rather than the 
rock according to simulations. This effect has been accounted for in the expected flux calculations. 

Stopping muons with track length > 7 m (~1.6 GeV) in the ID are selected for further analysis. The 
stopping muon track length is determined by calculating the muon momentum from a photoelectron count 
in the same way as in the contained event analysis ||3|]. This cut eliminates short tracks that are very 
close to the PMTs and thus are hard to reconstruct, provides an energy threshold, and also eliminates 
nearly all contamination of the upward stopping muon signal by pions. These pions are photoproduced 
by cosmic-ray muons outside the detector and are observed to have a soft spectrum ^ , such that residual 
contamination upward-going pions meeting the 7 m track length requirement is estimated to be < 1%. 
The nominal detector effective area for upward-going muons with a track length > 7 m in the ID is 
~1200 m2. 
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137 upward stopping muons which satisfy a cos O < cut are found, where Q is the zenith angle of 
the muon track, with cos0 = — 1 corresponding to vertically upward-going events. 

Details of the muon track reconstruction method and data reduction algorithm are similar to those 
of through-going muons described elsewhere [|]. The total detection efficiency for the complete data 
reduction process for upward stopping muons is estimated by a Monte Carlo simulation to be >98% 
over — 1 < COS0 < 0. The validity of this Monte Carlo estimate is in turn checked with real cosmic-ray 
downward stopping muons, taking advantage of the up/down symmetry of the detector geometry. 

Because of multiple Coulomb scattering and finite muon fitter angular resolution (~ 1° [Q), some of the 
abundant downward-going cosmic-ray stopping muon background may be reconstructed with cos O < 
and contaminate the neutrino induced upward stopping muon sample. Figure ^ illustrates the estima- 
tion of this contamination by extrapolation of the downward-going zenith angle distribution to the flat 
neutrino- induced signal near the horizon. This background falls exponentially with decreasing cos0, 
the contribution to apparent upward stopping muons is estimated to be 10.8±4.2 events in the zenith 
angle bin with cos0 between -0.2 and 0. This is a larger contamination than was present in the upward 
through-going muon sample due to the lower muon energies (and correspondingly larger multiple scat- 
tering) of the stopping muon sample, and is more uncertain due to lower statistics. As an independent 
crosscheck of the background contamination, the energy spectrum of the observed stopping muons is 
shown in Fig. ^. That shape is consistent with the expected and no significant contamination is observed 
at the 1.6 GeV muon energy analysis threshold. 

To analytically calculate the expected upward stopping (and through-going as well) muon flux, em- 
ployed is the combination of the Bartol atmospheric neutrino flux model jw^ , a neutrino interaction model 
composed of quasi-elastic (QE) scattering + single-pion (Itt) production ||l^ -I- deep inelastic (DIS) 
scattering multi-pion production based on the parton distribution functions (PDF) of GRV94 |]l^ with 
the additional kinematic constraint of Vl^ > 1.4 GeV/c^ (where W is the invariant mass of the hadronic 
recoil system), and Lohmann's muon energy loss formula in standard rock jl^. 

This expected flux is compared to three other analytic calculations to estimate the model-dependent 
uncertainties of the expected muon flux. The other flux calculations use various pairwise permutations 
I or the atmospheric neutrino flux model calculated by the Bartol group |l^, the 
GRV94DIS PDF or the CTEQ3M ||l|lPDF. This comparison yields ±10% of difference in the overaU flux 
and -2% to -1-1% for the bin-by-bin shape difference in the zenith-angle distribution. The shape difference 
is due mostly to differences in the input flux models. 

The expected muon flux f&j^gQ resulting from the above calculation is 0.73 ± 0.16 x lO^^'^cm^-^s^^sr"^ 
(cos 8 < 0), where the estimated theoretical uncertainties are described in Table || The dominant 
error comes from the overall normalization uncertainty in the neutrino flux, which is estimated to be 
approximately ±20% Q^lljl^] above several GeV. 

Given the detector live time T, the effective area for upward stopping muons S{Q), and the detection 
efficiency e(6), the upward stopping muon flux is calculated using: 

^1 11 

^^siQ,) S{e,)2n T 

where the index j runs over observed events, 27r is the total solid angle covered by the detector for upward 
stopping muons, and N is the total number of observed muon events (137). Subsequently, we subtract the 
cosmic-ray muon contamination (10.8 events) from the most horizontal bin (-0.2<cos6<0). The resulting 
observed upward stopping muon flux is: = 0.39 ± 0.04(stat.) ± 0.02(sys.) x 10~^^cm~^s~^sr~^. If 

instead of subtracting this background, we simply omit upward stopping muons coming from the thin-rock 
region [Q| in the most horizontal cos0 bin (-O.l<cos0<O and 60° < cf) <310°), the calculated flux differs 
from the background-subtracted flux by -1.1%. Systematic errors in the experimental measurement are 
summarized in Table 

The flux as a function of zenith angle, {d^^*°P / dfl) , is shown in Fig. |^. Due to limited statistics for 
stopping muons, 5 angular bins are used instead of the 10 bins used for Super-K through-going muons. 
With the present statistics, the shape of the normalized distribution is consistent with the no-oscillations 
hypothesis (x^/d.o.f. = 4.1/4 corresponding to 39% probability). However, the overall flux of upward 
stopping muons observed is substantially depressed from that expected, as this shape comparison is made 



of the Honda flux |15 
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after multiplying the expected flux by a free-running overall normalization factor (l+a) whose best fit 
value is a = —51% calculated by a shape fit . This can be compared with an observed no-oscillations 
a = —14% in the through-going muon sarnple 0]. 

Adding the upward through- going flux |^] <I)*"^"^ we obtain a detector-independent total upward-going 
muon flux of 2.13±0.08(stat.)±0.03(sys.) x lO-^^cm-^s-isr"! with muon energy > 1.6 GeV. 

The zenith angle distribution of the summed flux, (c?<I>''*°P+*'"'"/c?ri), is shown in Fig. |[ 

A more useful physical quantity than the absolute flux for probing neutrino oscillations is the 
stopping/through-going flux ratio H = ^"^ov j ^i^ru ^ which cancels much of the large (~ 20%) uncer- 
tainty in the neutrino flux normalization and the neutrino interaction cross sections |18| . The systematic 
theoretical and experimental uncertainties in "R, are summarized in Table [I] and Table IV, respectively. 
The measured 7^ is 0.22 ± 0.02(stat.) ± 0.01(sys.), while the expected Tltheo is 0.37t^(theo.). The 
probability that the observed ratio could fluctuate this far below the expectation is 0.39%. The zenith 
angle distribution of the ratio is shown in Fig. ^ 

Using these data, we derive probability contours on the neutrino oscillation parameter {sin^20, Am^) 
plane for the i'^ <-*■ oscillation hypothesis, as shown in Fig. |^, based on a defined by: 



Y = mm 

/3 



d{stop) 




W - ni^^^{sin^20, Ato2) X (1 + /?) 



where TV is the observed ratio in the z-th cos0 bin, CTt^.^ ^ the experimental statistical error, a-j^i 
the bin-by-bin uncorrelated systematic error (~2%) estimated by adding uncorrelated theoretical and 
experimental systematic errors in Table III and Table IV in quadrature, and T^J/jg^ the expected ratio. 
The P for 'R-lf^,,,, represents a stopping-to-throughgoing relative muon flux normalization factor with error 
(7^=14% estimated by the correlated theoretical and experimental systematic errors in Table [II and 
Table IV added in quadrature. 



The allowed region thus obtained in Fig. R is in good agreement with that found in the Super-K 
contained event analysis 1^ . The minimum x locations on the Am^ — sin^ 26 plane and corresponding 
parameter values for various conditions are listed in Table As the minimum lies in the unphysical 
region, the contour is drawn according to the prescription for bounded physical regions given in Ref . . 
Because the surface has a rather broad minimum, the specific best-fit oscillation parameter values 
cited are of less importance than the contours shown in Fig. |^. If we replace the Bartol neutrino fiux |0 
by the Honda's fiux Q and/or the GRV94DIS parton distribution functions ^ by CTEQ3M ||l6|, the 
allowed region contours do not change significantly. 

Finally, we performed an oscillation analysis which use all the upward-going muon information available 
by simultaneously fitting to the upward through-going and stopping muon zenith angle distributions. In 
this analysis, is defined by the sum over 10 upward-throughgoing and 5 upward-stopping muon zenith 
angle bins: 



Y — mm 



-1 V ^i'^ltat)' + i<ysr 



E 



where is the observed muon flux in the i-th cosO bin, al^^f the experimental statistical error, cr^j,,, 

(~2 to ~4%) the bin-by-bin uncorrelated theoretical and experimental systematic errors in Table [| and 
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Table ^ added in quadrature for stop or from Ref. for thru , the expected muon flux, a the 

running overaU flux normahzation factor with error CTq of 22% estimated by adding in quadrature the 
correlated theoretical systematic errors in Table | /3 the running stopping-to-throughgoing relative muon 
flux normalization factor with error (T^=14% estimated by the correlated theoretical and experimental 
systematic errors in Table HI and Table IV added in quadrature. The results are shown in Fig. |^ and are 
consistent with those from the Super-K contained event analysis. As is shown in Table 0, the minimum 
falls in the unphysical region so the contours are drawn according to the prescription for bounded 
physical regions given in Ref. |19|. 

The contamination due to v,, charged current interactions and neutral current interactions in the 
upward-going muons is estimated to be < 1% by a Monte Carlo simulation and is neglected in these 
analyses. The contribution of possible Vr interactions in the rock below to the upward stopping muon 
flux is suppressed by branching ratios and kinematics to less than a few percent, and is also neglected. 

The observed overall upward stopping muon flux alone and its zenith angle distribution do not conflict 
with the expected no-oscillations values signiflcantly within the present statistical and systematic errors. 
However, in order to simultaneously explain both the stopping and through-going upward muon fluxes 
using both an analysis based on the stop/through ratio TZ and a combined fit of upward stopping and 
through-going muon zenith angle data, the no-oscillations theoretical expectations reproduce the observed 
fluxes poorly at the 0.39% C.L. (stop/thru ratio averaged over zenith angle) and 0.87% C.L, respectively, 
while the i^^ <-> i't oscillation assumption is consistent with observations. This result supports the evidence 
for neutrino oscillations given by the analysis of the contained and through-going muon atmospheric 
neutrino events by Super-K. 

We gratefully acknowledge the cooperation of the Kamioka Mining and Smelting Company. The Super- 
Kamiokande experiment has been built and operated from funding by the Japanese Ministry of Education, 
Science, Sports and Culture, and the United States Department of Energy. 



TABLE I. List of theoretical uncertainties in upward stopping muon flux calculation. 



Error 


source 


Error 


typical 1/ flux normalization 


±20%=^ 


Choice of ly flux/PDF 








overall flux 


±10%" 




bin by bin 


-2% to +1% 



"Theoretical cosB bin-by-bin correlated uncertainty 
''Theoretical cosO bin-by-bin uncorrelated uncertainty 



TABLE IL List of experimental systematic errors in upward stopping muon flux measurement. 



Error source 


Error 


7 m track length cut 


±5%" 


Live time 


±1%" 


Reconstruction efliciency 





"Experimental cosB bin-by-bin correlated systematic error 
''Experimental cosO bin-by-bin uncorrelated systematic error 

Note that the uncertainty of the background subtraction is included in the statistical error. 
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TABLE III. List of theoretical uncertainties in stopping/through-going muon flux ratio. 



Error source 


Error 


Primary spectral index ±0.05 


±13%" 


Choice of v flux 




overall ratio 


±1%" 


bin by bin 


-1% to ±2%'= 


V cross section 


± 4%" 



"Theoretical cos0 bin-by-bin correlated uncertainty 
''Theoretical cosO bin-by-bin uncorrelated uncertainty 

'^QE, Itt and DIS cross sections are changed independently by ±15%, and the effect of this on Tl is noted here. 



TABLE IV. List of experimental systematic errors in stopping/through-going muon flux ratio. 



Error source 


Error 


7 m track length cut 


±5% " 


Live time 


±1% " 


stop/thru misidentification 


~ ±1% " 


Reconstruction efficiency 


~ ±1% 



"Experimental cosO-bin-by-bin correlated error 
''Experimental cosO-bin-by-bin uncorrelated error 



TABLE V. Summary of fit results. See text for definitions of a and /3. 



Case 




single 


a 


/3 


xi„iv/DOF 


Probability 


Flux vs. zenith angle (Fig. 5): 






-0.51 




4.1/4 


39% 


"R. analysis (Fig. ^): 














Fit in the physical region 


3.1 X 10-^eV^ 


1.00 




-0.084 


0.76/3 


86% 


Fit over all space 


3.9 X 10-^eV^ 


1.19 




-0.003 


0.30/3 


96% 


No oscillations 








-0.35 


8.4/5 


13% 


(stop±through) cos(G) analysis (Fig. ^): 
Fit in the physical region 


3.9 X 10"^eV^ 


1.00 


0.061 


-0.083 


8.8/13 


79% 


Fit over all space 


3.9 X 10"^eV^ 


1.13 


0.10 


-0.012 


8.2/13 


83% 


No oscillations 






-0.16 


-0.36 


31.0/15 


0.87% 
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FIG. 1. The energy spectra of the parent neutrinos of upward stopping muons (left) and upward through-going 
muons (right). The dashed lines are the result of using the Bartol input fluxes, and the solid lines are from 
Honda's calculations. 



8 



10 



(|)<60°,(t)>310° 



o 60°<(^<310° 



S 10 



= 10 r 




1 r 



COS© 



FIG. 2. Zenith angle distribution of stopping muons near the horizon observed by Super-K. Filled triangles 

(open circles) indicate events coming from direction where the rock overburden is thick > 15000 m.w.e. (shallow 
> 5000 m.w.e.). The two distributions are normalized to a common azimuth angle range (0 - 360 degrees). 
The solid line is the zenith fit to the shallow rock events used to estimate the cosmic-ray muon background 
contamination in the —0.2 < cos© < bin. 
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FIG. 3. The energy distribution of upward stopping muons, illustrating the freedom of the data both from 

threshold effects due to the 7 m (1.6 GeV) muon pathlength cut and from upward-scattered low energy stopping 
cosmic-ray muons. The solid histogram is the normalized distribution expected from the no-oscillations hypothesis. 
The normalization is made so that the integrated number of expected events with muon energy > 1.6 GeV (the 
vertical line) corresponds to the observed number of events (=137). Preselection is applied at muon energy 
~1.4 GeV. 
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^ 4 




COS0 

FIG. 4. Upward stopping (filled circles) and through-going (open circles) muon fluxes observed in Super-K as 
a function of zenith angle. The error bars indicate uncorrelated experimental systematic plus statistical errors 
added in quadrature. The lower (upper) solid histograms show the expected upward stopping (through-going) 
muon flux based on the Bartol neutrino flux without oscillation. Also shown as lower (upper) dashed histograms 
are the expected stopping (through-going) muon flux assuming the best flt parameters of the combined analysis in 
the physical region at (sin^ 261, Am^) = (1.0, 3.9 x 10"^eV^), a = -1-0.061 and /3=-0.083 for the ^ Vr oscillation 
case. 
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FIG. 5. Upward stopping+through-going muon flux (fiUed circles) observed in Super-K as a function of zenith 
angle. The error bars indicate uncorrelated experimental systematic plus statistical errors added in quadrature. 
The solid histograms show the expected flux based on the Bartol neutrino flux without oscillation. Also shown as 
dashed histograms are the expected flux assuming the best flt parameters of the combined analysis in the physical 
region at (sin^ 26, Am^) = (1.0, 3.9 x 10"^eV^), a = +0.061 and ^=-0.083 for the v^, *-* Vr oscillation case. 
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FIG. 6. Upward stopping/through-going muon flux ratio (filled circles) observed in Super-K as a function 
of zenith angle. The error bars indicate uncorrelated experimental systematic plus statistical errors added in 
quadrature. The solid histograms show the expected the ratio based on the Bartol neutrino flux without oscillation. 
Also shown as dashed histograms are the expected ratio assuming the best flt parameters of the combined analysis 
in the physical region at (sin^ 20, Am^) = (1.0,3.9 x 10"^eV^), a = +0.061 and /3=-0.083 for the ^ v-r 
oscillation case. 
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FIG. 7. The allowed region contours at 68% (dotted contour), 90% (thick solid), and 99% (dashed) C.L. ob- 
tained by the combined analysis of Supor-K upward stopping and through-going muons drawn on the (sin'^ 29, Am^) 
plane for ^ Vr oscillations. The star indicates the best fit point at (sin'^ 29, Am^) — (1.0, 3.9 x 10"'^ eV^) in the 
physical region. The allowed region contour indicated by solid thick labelled line with "STOP/THRU" is made 
based on the Super-K stopping/through-going muon ratio alone at 90% C.L. Also shown is the allowed region 
contour (the remaining solid thin line) at 90% C.L. by the Super-K contained event analysis. The allowed regions 
are to the right of the contours. 
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